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Abstract 

Any unitary transformation of quantum computational networks is explic- 
itly decomposed, in an exact and unified form, into a sequence of a limited 
number of one-qubit quantum gates and the two-qubit diagonal gates that 
have diagonal unitary representation in usual computational basis. This de- 
composition may be simplified greatly with the help of the properties of the 
finite-dimensional multiple-quantum operator algebra spaces of a quantum 
system and the specific properties of a given quantum algorithm. As elemen- 
tary building blocks of quantum computation, the two-qubit diagonal gates 
and one-qubit gates may be constructed physically with one- and two-body 
interactions in a two-state quantum system and hence could be conveniently 
realized experimentally. The present work will be helpful for implementing 
generally quantum computations with any qubits in those feasible quantum 
systems and determining conveniently the time evolution of these systems in 
course of quantum computation. 



1. Introduction 

Since it has been discovered that quantum computers can be much more 
powerful than their classical counterparts (Feynman 1982; Deutsch 1985; 
Shor 1994), it becomes of great practical importance to realize the quan- 
tum computers. A variety of quantum systems have been explored to build 
such quantum computers such as trapped ions (Cirac & Zoller 1995), nu- 
clear spins in molecules (Gershenfeld & Chuang 1997; Cory et al. 1997) and 
in solid states (kane 1998), and Josephson junction arrays in superconduc- 
tors (Makhlin et al. 1999), etc. Very recently, nuclear magnetic resonance 
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(NMR) (Ernst et al. 1987; Freeman 1997) is used to realize experimen- 
tally the Deutsch-Jozsa s algorithm (Deutsch & Jozsa 1992, Chuang et al. 
1998a) and the Grover s algorithm (Grover 1997, Chuang et al. 1998b, Jones 
et al. 1998). It is expected that in the short term quantum computation 
with as many as ten qubits will be implemented, although there are many 
problems that need to be solved such as how to overcome the effects of de- 
coherence and dephase in quantum systems and how to realize fault-tolerant 
and error-correlating quantum computation. Quantum computation should 
be performed within the characteristic time of decoherence and dephase in a 
quantum system. In practice, this requires that elementary building blocks 
of quantum computation should be chosen suitably. For example, theoreti- 
cally any quantum computation should be built exactly out of a sequence of 
the building blocks with a length as short as possible, while these building 
blocks can be exactly constructed theoretically and could be physically real- 
ized conveniently in a feasible quantum system. Quantum gates were firstly 
suggested by Deutsch (1989) as the elementary building blocks to construct 
any quantum computational networks. A universal quantum gate, by copying 
itself and then wiring together, suffices to construct any unitary transforma- 
tion of quantum computation. It has been shown theoretically (Deutsch et 
al. 1995, Lloyd 1995) that almost every quantum gate that operates on two 
or more qubits is universal, but the universal quantum gates that are really 
considered as conveniently realizable gates are the three-qubit gates, e.g., 
Toffolis (1981), Fredkins (1982), and more general Deutsch s (1989) gate 
and the simpler two-qubit gates (Barenco 1995 & Sleator and Weinfurter 
1995). Quantum gates with more qubits may not be universally attractive 
since their construction and implementation may be usually complicated in 
quantum systems and the theoretical construction of any unitary transfor- 
mation of quantum computation out of these gates is usually not carried 
out easily. This may be the main reason why many investigators (DiVin- 
cenzo 1995a; Barenco 1995 & Sleator and Weinfurter 1995) have suggested 
the simpler universal quantum gates with only two qubits as the building 
blocks. Barenco et al. showed further that a set of quantum gates that con- 
sists of all one-qubit gates and the two-qubit XOR gate suffices to build any 
unitary transformation, although one-qubit gates and the simpler two-qubit 
XOR gate are not universal. However, the present theoretical composition 
of the three- or two-qubit universal gates or even the XOR gate along with 
one-qubit gates to form quantum networks usually disregards some useful 
properties of quantum systems and the specific properties of a given quan- 
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turn computation. As a consequence, it needs probably an infinite number of 
such quantum gates to build exactly a given unitary transformation of quan- 
tum computation. Obviously, this is impractical for quantum computation 
to be performed in a quantum system in which the effects of decoherence and 
dephase are not negligible. 

In this paper any unitary transformation of quantum computational net- 
works is decomposed explicitly as a sequence of a limited number of one- and 
two-body elementary propagators, i.e., one-qubit gates and the two-qubit 
diagonal gates Rki(^ki) = Gxp(—iXki2Ik z Ii z ) that have diagonal unitary rep- 
resentation in the conventional computational basis, in an exact and unified 
form. Quantum computation with any qubits is then implemented by per- 
forming the unitary transformation of a sequence of these one- and two-qubit 
gates on the input quantum state. In constrast to the usual composition of 
quantum gates to form quantum networks here is emphasized on the de- 
composition of a given quantum network as a sequence of the elementary 
building blocks. This decomposition may be simplified greatly with the aid 
of the properties of the finite-dimensional multiple-quantum operator alge- 
bra spaces of quantum systems (Miao, 2000a) and the specific properties of 
a given quantum algorithm. That the two-qubit diagonal gate Rki{^ki) is 
chosen as an elementary building block stems from several considerations. 
Firstly, such choice for building blocks is beneficial to the exploitation of 
the properties of the finite-dimensional multiple-quantum operator algebra 
spaces of quantum systems to simplify the decomposition. Secondly, the 
two-qubit diagonal gates supplemented with all one-qubit gates suffice to con- 
struct exactly any unitary transformation of quantum computation. Thirdly, 
from the physically realizable point of view one-qubit gates should be the sim- 
plest gates, while according to matrix (operator) algebra properties quantum 
gates that have diagonal unitary representation matrices in the conventional 
computational bases {|00...00), |00...01), |ll...ll)} are the elementary and 
simple gates. Particularly, the two-qubit diagonal gate Rm(Xm) that has 
diagonal unitary representation matrix: 



should be the most elementary and simplest building blocks, although it 
is not a universal gate. Moreover, the two-qubit diagonal gate Rki(^ki) can 
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be conveniently built up with one- and two-body interactions such as the gen- 
eral neighbor interaction and could be easily realized in a quantum system 
with N two-state particles such as coupled multispin systems in molecules or 
in solid states, trapped ions, superconducting Josephson junction arrays, etc. 
Fourthly, the two-, three-, and N-qubit (N>3) universal gates and even the 
XOR gate can be expressed exactly as a simple sequence of the two-qubit 
diagonal gates along with one-qubit gates. 

2. The decomposition of unitary transformation 

Quantum computation is a reversible process (Bennett, 1973). It can be 
thought of as a unitary transformation acted on the input state and obeys 
the laws of quantum mechanics (Benioff, 1980; Deutsch, 1985; DiVincenzo, 
1995b). The time evolution of a quantum system from the initial state to 
the output state during quantum computation then can be described by a 
time-evolutional propagator U (t) that obeys the Schrodinger equation: 

jU(t) = -iH(t)U(t) (h=l) (1) 

where H(t) is the effective Hamiltonian for the system to perform the quan- 
tum computation. The effective Hamiltonian H(t) characterizes generally the 
specific properties of the quantum computation. In general, quantum com- 
putational networks are composed of a sequence of quantum circuit units 
(Deutsch, 1989). Each such circuit unit performs the unitary transforma- 
tion with a propagator U k {t k ) associated with the time-independent effective 
Hamiltonian H k in the interval t k , while the total propagator U(t) can be ex- 
pressed as a sequence of the propagators U k (t k ). Then it follows from Eq.(l) 
that 

U(t) =n U k (t k ) =n exp(-iH k t k ) (t =£ t k ) (2) 

k k k 

The quantum computation then can be implemented by acting a sequence of 
the propagators U k (t k ) on the input state. Therefore, it becomes clear that 
the problem to be solved is how to exactly decompose theoretically the prop- 
agators U k (t k ) of quantum circuit units as a sequence of a limited number of 
one- and two-qubit gates and how to build up these simple gates experimen- 
tally in an accessible quantum system. Not loss of generality, the quantum 
system is considered as a physical system consisting of N two-state particles. 
This system may be nuclear spins in molecules or in solid state, trapped 
ions, and superconducting Josephson junctions, etc. Here for simplification 
the complete decomposition of the propagators is described explicitly in a 
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coupled spin (1=1/2) system, which is formed by N two-state nuclei with 
magnetic quantum number 1=1/2. 

The effective Hamiltonian Hj associated with each quantum circuit unit 
can be generally expanded as a linear combination of base operators {-£>&} of 
the Liouville operator space of the spin system (Ernst, et al. 1987): 

Hj —J2 a k B k (3) 

k 

As suggested recently (Miao, et al. 1993 & 1997; Miao, 2000a), to determine 
exactly and analytically time evolution of the spin system the propagator cor- 
responding to this Hamiltonian is first decomposed into an ordered product 
of a series of elementary propagators 

Uj{tj) = expi-iHjtj) =n Rs(Xs) (4) 

s 

The elementary propagator is defined by 

R S (K) = exp(-i\ s B s ) (5) 
where A s is a real parameter and B s a Hermite base operator. Obviously, 
the elementary propagator is also a quantum gate. Actually, the decom- 
position of Eq.(4) can be achieved in an exact and unified form. Firstly, 
the propagator Uj(tj) is converted unitarily into a diagonal unitary opera- 
tor, which has diagonal unitary representation in usual computational basis, 
by making a sequence of elementary unitary transformations. Then each 
such elementary unitary transformation and the diagonal unitary operator 
are further decomposed into a product of a series of elementary propagators, 
respectively. The decomposition of Eq. (4) can be further simplified with the 
help of the properties of the Liouville operator spaces and its three subspaces 
(Miao, 2000a): the even-order multiple-quantum, the zero-quantum, and the 
longitudinal magnetization and spin order operator subspace. When the ef- 
fective Hamiltonian Hj is a member of the longitudinal magnetization and 
spin order operator subspace, the propagator Uj(tj) is simply expressed as 
a sequence of elementary propagators built up with the base operators of 
the subspace (see below). If Hj is a member of the zero-quantum operator 
subspace, one first makes a zero-quantum unitary transformation on Uj(tj) 
to convert it into the diagonal unitary operator and then further decom- 
poses the zero-quantum unitary operator and the diagonal unitary operator 
as a sequence of elementary propagators, respectively. When the effective 
Hamiltonian Hj is a member of the even-order multiple-quantum operator 
subspace, one makes the even-order multiple-quantum and subsequently the 
zero-quantum unitary transformation on Uj(tj) to convert it into the diago- 
nal unitary operator. The even-order multiple-quantum, the zero-quantum, 
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and the diagonal unitary operator can be further decomposed as a sequence 
of elementary propagators, respectively. If Hj is not a member of any one of 
the above three subspaces but a member of the Liouville operator space, one 
first converts it unitarily into a member of the even-order multiple-quantum 
operator subspace by making an odd-order multiple-quantum unitary trans- 
formation on the Hamiltonian, then a further decomposition for the propaga- 
tor Uj(tj) can be carried out with the help of the properties of the even-order 
multiple-quantum operator subspace. 

It is clearly shown from the closed property of operator algebra space 
that any quantum gate built up with an arbitrary operator of any one of the 
three aforementioned operator subspaces is a non-universal gate. These gates 
can form another set of non-universal gates that may be different from one- 
qubit gates and collection of one-qubit gates and the classical gates (Deutsch, 
1995). 

On the basis of the decomposition of Eq.(4) time evolution of a system in 
the course of quantum computation can be determined directly by acting the 
decomposed propagator on the input state in a quantum system or on the 
initial density operator in a quantum ensemble with the help of the rotation 
transformation between any two base operators. The rotation transformation 
can be generally derived from the Baker-Campbell-Hausdoff formula (Ernst, 
et al. 1987): 

R^BMK)- 1 = £ { -^C n (6) 

n=0 

where Co = B r and C n = [B s , C n _i] and particularly, if [B s , [B s , B r }} = aB r 
the transforamtion (6) reduces to a simpler closed form 

R s (X s )B r R s (X s )- 1 = B r cos(v^A s ) - ^[B s , B r \ sin(^A s ) (7) 
For a coupled N-spin (1=1/2) system the proper base operators {B^} 
of the Liouville operator space are usually chosen as the Cartesian product 
operators (S^rensen, et al. 1983; Ernst, et al. 1987): 

{Bk} = {E, h ia , 2Ikiah 2 p, — ■> 2 n ~ 1 J fclQ J fc2/3 ...J fcn< 5; 1 < n < N} (8) 
where E is unit operator and Ik ia («, (3,8 — x, y, z) are spin angular momen- 
tum operators for the kith, spin in the system (J^ = \ok^ ® is the Pauli s op- 
erator). Such direct product operator set contains any n-body (N > n > 1) 
interaction terms {2 n ~ 1 I ka Ii/ 3 ...I m s}. It follows from Eqs.(2), (4), and (5) 
that the propagator Uj(tj) is usually expressed as a sequence of elementary 
propagators built up with any n-body (N > n > 1) product operators in set 
(8). Actually, the propagator Uj(tj) can be further expressed as a sequence of 
the elementary propagators built up only with one- and two-body operators 
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of the set (8). 

Each base operator of set (8) can be converted unitarily into a member 
of the longitudinal magnetization and spin order operator subspace (Miao, 
2000a), where the base operators of the subspace are usually chosen as the 
longitudinal magnetization and spin order product operators in the N-spin 
(1=1/2) system (Miao, et al. 1993; Miao, 2000a): 

z Imz i • • • j 2 

zhz---^Nz } (9) 

Then any elementary propagator defined by Eq.(5) can be converted unitar- 
ily into an elementary propagator built up with a product operator of the 
subspace by applying a limited number of 90 degree electromagnetic pulses. 
A typical example is shown below 

exp(zf ifcy) exp(-if ite) exp(-ifl mx ) 

J-kx-l-ly-lmy ^kz^lz---^mz 

where the unitary transformation B = UAU + is denoted briefly as A— — B. 

Therefore, the basic building blocks for the propagator Uj(tj) are those ele- 
mentary propagators built up with the base operators of the subspace and 
the one-body elementary propagators of Eq.(5). On the other hand, any 
elementary propagator constructed with an n-body (N > n > 1) product 
operator of the subspace can be readily decomposed as a product of a series 
of the elementary propagators built up only with the two-body product op- 
erators in set (9) and the one-body base operators in set (8). This can be 
achieved by utilizing recurrently the following decomposition: 
exp(-i\2 n I klZ ...I knZ I kn+lZ ) = 

V n e W (-i\2 n - 1 I klZ ...I kn _ lZ I kn+lZ )V+ (n > 2) (10) 

where 

K = exp(-qi kn+ix ) exp(-inI knZ I kn+lZ ) exp(if I kn+lX ) exp(-zf I kn+iy ) 
As a consequence, it follows from Eqs.(4) and (5) that the propagators Uj(tj) 
of the quantum circuit units and hence the total propagator U(t) of quan- 
tum computation can be decomposed completely into a product of a series 
of one-body elementary propagators and the two-body diagonal elementary 
propagators R k i(X k i) built up with the product operators {2I kz Ii z }. 

Evidently, any operator of the longitudinal magnetization and spin order 
operator subspace has the diagonal representation in usual computational ba- 
sis and any two base operators B™ and B™ of the subspace are commutable 
with each other. Then any diagonal operator of the system, i.e., an oper- 
ator that has diagonal representation in usual computational basis, can be 
expressed as a sum of the base operators of the subspace (Miao, 2000a). If 
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the effective Hamiltonian Hj (=^ Ofc-B™) associated with a quantum circuit 

k 

unit is a member of the subspace, the corresponding propagator exp(-iHjtj) 

can be readily decomposed as a product of a series of elementary propagators 

constructed with the base operators of the subspace: 

exp(-iHjtj) =n exp(-ia k B^tj) (11) 
k 

Equation (11) is very useful for the decomposition of the total propagator of 
a given quantum algorithm. 

3. Preparation of the elementary building blocks 

The two-qubit diagonal quantum gate R k i(X k i) could be easily prepared 
in many two-state physical systems. As an example, its preparation is de- 
scribed explicitly in an accessible coupled N-spin (1=1/2) system. In general, 
the external electromagnetic field such as radiofrequency (RF) field is used to 
control the process of quantum computation and in the coupled spin (1=1/2) 
system with ravelled resonances each one-body elementary propagator de- 
fined by Eq.(5) may be prepared by utilizing selective pulses (the weak RF 
field) (Freeman, 1997). The spin Hamiltonian for the system in a strong 
static magnetic field is written as (Ernst, et al. 1987) 

H =£ n k I kz + £ nJ kl 2I kz I lz (12) 

k k<l 

where it is assumed that the internuclear interaction is weak with respect 
to the Zeeman interaction and the interaction between the system and its 
environment, which results in decoherence and dephase, is negligible. This 
Hamiltonian that consists of one-body {I kz } and two-body {2I kz Ii z } inter- 
actions is responsible for preparing experimentally the two-qubit diagonal 
quantum gates R k i(X k i). Figure 1 presents the quantum circuit unit (the 
NMR pulse sequence) for the preparation of the elementary propagator built 
up with the direct two-body interaction between two spins k and I, where 
spin echoes refocus all the undesired one- and two-body interactions and only 
leave selectively the desired two-body interaction 21 kz I\ z in the Hamiltonian 
(12) by combining selective 180 degree pulses. If there is not direct inter- 
action between any two spins k and m, their indirect two-body interaction 
2I kz I mz may be achieved through a directly neighbor coupling network such 
as k — I — ... — s — t — m in the system: 

exp(-inl kx l rx ) exp(-inl ky l ry ) _ 

^^kz^mz ^i-rz^-mz ^^sz^-rnz 
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exp(-iTrI sx I tx ) exp(-iirl sy l ty ) 

^^tz^mz 

For quantum dots (Barenco, et al. 1995a; Loss & DiVincenzo, 1998) the 
diagonal gates Rki(^ki) could be prepared in an analogous way to the above 
approach. In trapped ion system (Cirac & Zoller, 1995) the diagonal gates 
could be implemented by six laser pulses (see Appendix B) and in supercon- 
ducting Josephson junction arrays (Makhlin, et al. 1999) they could also be 
prepared easily (Miao, 2000b). 

4. Application to the universal quantum gates and quantum al- 
gorithms 

It is easy to carry out the explicit decomposition of the total propagators 
for N-qubit quantum algorithms such as the Deutsch-Jozsa, Grover, quantum 
Fourier transform algorithm, etc. and for the two-, three-, and N-qubit (N>3) 
universal quantum gates. Several typical examples are given explicitly below. 

4.1 The two-, three-, and any N-qubit universal quantum gates 

The unitary representation matrix Un of the N-qubit universal gate ( 
Deutsch, 1989; Barenco, 1995 & Barenco, et al. 1995b) can be generally 
written as 

U N = E+Diag{0, 0, 0, 1) (g)( 

where the matrix [uij] acting on the Nth. qubit is u{2) unitary matrix and 
can be generally expressed as 

U{2) = T N expH(^ + (p 1 I Nz )]T+ (14) 
where T/v — exp(— ial^ z ) exp(— i(3I^ y ). Then the unitary operation Un can 
be decomposed completely as a simple sequence of one-body elementary prop- 
agators and the two-body diagonal^ elementary propagators: 

U N = T N exp(-iH N t)T+ (t = 1) (15) 
where the diagonal operator H N is a member of the longitudinal magnetiza- 
tion and spin order operator subspace: 

N N N 

H N = ^0+ E ti k I kz + E J'kl^kzllz+ E J' klm ±h z llzlmz+-, (16) 
k=l l>k=l m>l>k=l 

here unit operator E is omitted. All the parameters in Eqs.(14)-(16) are de- 
termined directly from the elements {uij} of the matrix u(2) (Miao, 2000a). 

In particular, for the two- and three-qubit gates (Barenco, 1995 & Deutsch, 
1989) the diagonal unitary operators can be respectively written as 
exp(— iH 2 t) = exp(— iQ ) exp(— iQ^I^) exp(— iQ' 2 l2 Z ) 



+ 



un Ui2 

U 2 1 U 22 



(13) 
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x exp(-if 12 2I lz I 2z ) (17) 

and 

exp(-iH 3 t) = exp(-ifio) II exp(-in' k I kz ) f] exp(-ij' kl 2I kz I lz ) 

k=l l>k=l 

x exp(-ij' 123 AI lz I 2z I 3z ) (18) 
By using Eq.(lO) the last three-body elementary propagator on the right- 
hand side of Eq. (18) can be further decomposed as a sequence of six one-body 
elementary propagators and three two-body diagonal elementary propagators 
and hence the Deutsch s three-qubit universal gate is exactly decomposed as 
a sequence of six two-qubit diagonal gates Rki(Xki) and thirteen one-qubit 
gates and one constant phase factor. 

4.2 The Grover s quantum search algorithm. 

In most quantum algorithms the first step is the creation of a superposi- 
tion. This may be achieved by applying Walsh- Hadamard transform on the 
groundstate in a quantum system. Any n-qubit Walsh-Hadamard transform 
is constructed by the direct product of n single-qubit M matrices (Grover, 
1997) and can be expressed as a sequence of one-body elementary propaga- 
tors and a constant phase factor: 

W = M 1 ®M 2 ®...®M n 

n n 

= exp(i?f) exp(-in Y, ifcx)exp(-i§ I ky ) 

k=l k=l 

In addition to the Walsh-Hadamard transform the basic unitary operations 
required by the Grover s algorithm (Grover, 1997) are the conditional phase 
shift operations represented respectively by the diagonal unitary matrix C 
and R: 

Cij = 0, if j; Cu = -1; if i = s; C« = 1, if s 
and Rij = 0, if i ^ j; Ru = -1, if i = 1; R i{ = -1, if i ^ 1 
The unitary operation C can be further expressed in the form of exponential 
operator 

C = exp(-iH c t) (t = 1) 
where the representation matrix elements of the diagonal operator H c can be 
derived from the matrix C as 

(H c )ij = 0, if i ^ j; (H c )ii = n, if i = s; (H c )u = 0, if i ^ s 
Therefore, the operator H c is a member of the longitudinal magnetization 
and spin order operator subspace and can be expressed as a sum of the base 
operators of the subspace. It takes the same form as the diagonal operator 
if at of Eq.(16) but with different parameters determined from the matrix 
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elements {{H c )ij}. As a result of Eq.(ll), the unitary operation C can be 
decomposed into a product of elementary propagators built up with the base 
operators of the subspace 

C = exp(-itto) exp(-iQ' k I kz ) Yl exp(-iJ M 2I kz Ii z ) 

k=l l>k=l 

x 11 exp(-iJ Hm 4I kz I lz I mz )... 

m>l>k=l 

According to Eq.(lO) the unitary operation C can be further decomposed 
completely as a sequence of one-body elementary propagators and the two- 
body diagonal elementary propagators. The diagonal phase rotation opera- 
tion R can be decomposed completely in an analogous way as the unitary 
operation C. Thus, each of the basic unitary operations W , C, and R (the dif- 
fusion transform D = WRW) in any n-qubit Grover s algorithm is expressed 
explicitly as a sequence of one-qubit gates and the two-qubit diagonal gates 
Rki{\i)- This result may be helpful to implement experimentally the algo- 
rithm with any qubits in an accessible two-state quantum system. 

4.3 The Deutsch- Jozsa s algorithm 

To decide certainly whether a function / : B n — ► B is balanced or 
constant, it needs to run unitary transformation Uf only once on the super- 
position (Deutsch & Jozsa, 1992; Jozsa, 1998; Cleve, et al. 1998): 

U f--^=E l*i>[4=(|0>-|l»] — 



^2™ Xi€B" V% 

-±= E (-i) /( ^)[4(|0)-|i>)] 



Therefore, Uf has a diagonal unitary representation and can be expressed in 
the form of exponential operator: 
U f = exp(-iH f t) (t = l) 
where the diagonal operator Hf has the representation matrix elements: 

Obviously, any N-qubit unitary operation U / can be readily decomposed com- 
pletely in a similar way to the unitary operation C in the Grover s algorithm. 



5. Discussion 

It is widely believed that the three-qubit universal gates are sufficient to 
build any quantum computation (Deutsch, 1989), but several investigators 
(DiVincenzo, 1995a; Barenco, 1995; Sleator and Weinfurter, 1995) showed 
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that any three-qubit universal gates can be further expressed as a sequence 
of the two-qubit universal gates and hence the latters are more basic units 
in quantum computation. In the paper it is shown that any three- and 
two-qubit universal gates as well as the XOR gates can be expressed as a 
simple sequence of one-qubit gates and the two-qubit diagonal gates. These 
simple gates can be constructed with natural one- and two-body interactions 
such as neighbor interaction and could be readily realized experimentally 
in a two-state quantum system. Therefore, the two-qubit diagonal gates 
should be also proper elementary building blocks to construct conveniently 
any quantum computation physically. 

The decomposition for any unitary transformation of the quantum net- 
work of a given quantum algorithm into a sequence of one-qubit gates and 
the two-qubit diagonal gates provides a good scheme for the quantum algo- 
rithm to be programmed on a quantum computer. The effective Hamiltonian 
of a given quantum algorithm characterizes generally the specific properties 
of the quantum algorithm and the operator algebra structure of the effective 
Hamiltonian may decide how the decomposition is implemented conveniently. 
Therefore, the decomposition may be achieved conveniently with the help of 
the operator algebra structure of the effective Hamiltonian and the proper- 
ties of the Liouville operator space and its three operator algebra subspaces. 
The explicit decomposition for any unitary transformation of quantum com- 
putational networks into a sequence of one-qubit gates and the two-qubit 
diagonal gates in an exact and unified form will be helpful for implementing 
generally any N-qubit quantum computation in feasible quantum systems 
and determining conveniently the time evolution of these systems in course 
of quantum computing. 

The effective Hamiltonian of a given quantum algorithm may also char- 
acterize generally the complexity of a quantum algorithm. Provided that 
the effective Hamiltonian consists of local interations of a quantum system 
subjected to the quantum algorithm, there is certainly a quantum computa- 
tional network that can simulate efficiently the quantum computation (Lloyd, 
1996). If a classical algorithm, which may not be efficient, is designed to solve 
an NP-problem in a classical digital computer and it can be translated into a 
quantum algorithm by replacing irreversible logic gates with the correspond- 
ing reversible gates according to the Bennett's suggestion (Bennett, 1973), 
now one wants to ask: can the NP-problem be solved efficiently with the 
quantum algorithm on a quantum computer? Evidently, this is impossible 
(Deutsch, 1985 & 1989). Is there other quantum algorithm to solve efficiently 
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the same NP-problem? If a quantum computational network is designed ac- 
cording to the mathematical structure and characteristic of the NP-problem 
and the quantum mechanical laws and if the effective Hamiltonian of the 
quantum network of the quantum algorithm is local, the network can solve 
efficiently the NP-problem. 
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Appendix A 

The multiple-quantum operator algebra spaces 

A p-quantum operator Q p is defined by (Miao, et al. 1993, Miao, 2000a) 

I z Q p \^ r ) = {M r + P)Q p \^ r ) ' {H = 1) (Al) 

where the wavefunction \ty r ) is an arbitrary eigenstate of the z-component I z 
of the total spin angular momentum operator of a spin system with its own 
eigenvalue M r 

I s \V r ) = M r |^ r ) (A2) 
The operator I z is also called the total magnetic quantum operator or the 
total longitudinal magnetization operator of the system. The definition of 
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Eq.(Al) of a p-quantum operator is general and independent of energy eigen- 
states of the system, although the wavefunction \ty r ) is also an eigenfunction 
of spin Hamiltonian of the system when the contribution of Zeeman interac- 
tion to the spin Hamiltonian is dominating. However, when spin Hamiltonian 
of a spin system contains non-secular interactions the operator I z usually 
does not commute with the spin Hamiltonian and in this case |\& r ) is not an 
eigenfunction of the spin Hamiltonian. 

The jo-quantum operator Q p has an explicit physical meaning that a new 
state Qp\^r) generated by Q p acting on an arbitrary eigenstate |\l/ r ) is also an 
eigenstate of the total magnetic quantum operator I z and its total magnetic 
quantum number raises p from the original one M r . It proves easily from the 
definition of the jo-quantum operator that the complete set of the p-quantum 
operators can construct a linear subspace of the Liouville operator space of 
the spin system since the sum of any two p-quantum operators is also a p- 
quantum operator. In particular, the complete set of zero-quantum operators 
is an operator algebra subspace of the Liouville operator space. This can be 
proven simply below. By expanding the eigenstates \^ r ) and Q P \^* T ) i n terms 
of the complete orthogonal and normalized eigenbase {\k)} of the operator 
I z with eigenvalues {M k }, respectively 

l*r) =E B rk (0)\k), for all k with M k = M r (A3) 
k 

and 

Qp\fr)=EBrk(0)C kl ( P )\l), (AA) 
k,l 

where sums run over all indexes I with Mi = (M r + P) and k with M k = M r , 
respectively, one can prove easily that the product operator of any two zero- 
quantum operators Qo a and is still a zero-quantum operator. It follows 
from Eq. (A4) that 

Q 0( sQ 0a \^r) = E B rk (0)C^0)C^(0)\m) (A5) 

k,l,m 

where sums run over all indexes k,l,m with M k , M;, M m = M r . Evidently, 
the product operator QopQoa is still a zero-quantum operator since all the 
eigenstates of the operator I z on the right-hand side of Eq.(A5) have the 
same eigenvalue equal to M r of \^ r ), 

Iz(QopQoa)\y r ) = M r (Q 0/3 Q 0a )\^ r ). 
Therefore, all the zero-quantum operators form an operator algebra sub- 
space of the Liouville operator space. As a direct result, the power oper- 
ator (Qo) n (n=l,2,...) of a zero-quantum operator Q Q is a zero-quantum 
operator and moreover, the exponential operator exp(±iAQo) of a Hermite 
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zero-quantum operator Q is a zero-quantum unitary operator and can be 
expressed as a sum of base operators {Qok} of the zero-quantum operator 
subspace 

exp(±UQ ) =E ^(Qo) n =E fk(±X)Qok (A6) 
k k 

where /^(zhA) are coefficients. There is an important property of the zero- 
quantum operator that any p-quantum operator does not change its quantum 
coherence order when it is acted on by an arbitrary zero-quantum operator. 
The proof for the property is simple. According to the definition Eq.(Al) of 
a p-quantum operator and Eqs.(A2)-(A4), one has 

QopQpQoal^r) = E B rk (0)C^(0)C lm {p)C^ n (0)\n) (A7) 

k,l,m,n 

where sums run over indexes k and I with = Mi = M r as well as m and 
n with M m = M n = (M r + P), respectively. Because all the eigenstates \n) 
on the right-hand side of Eq.(A7) have the same eigenvalue (M r + P), the 
state QopQpQoal^r) is an eigenstate of the operator I z and its own eigenvalue 
equals (M r + P), 

I z {QopQ P Qoa)\Vr) = (M r + P)(Q 0f3 Q p Q 0a )\y r ) . 
Therefore, the product operator QopQpQoa is a p-quantum operator, indi- 
cating that any p-quantum operator keeps its quantum coherence order un- 
changed when it is acted on by a zero-quantum operator. Particularly, any 
zero-quantum operator can be transferred into a sum of the base operators 
of the zero-quantum operator subspace by making a zero-quantum unitary 
transformation. This is really a direct consequence of the closed property of 
the zero-quantum operator subspace. 

In particular, it follows from the definition of Eq.(Al) of a zero-quantum 
operator that all the zero-quantum operators that are commutable with each 
other and also commute with the total magnetic quantum operator I z should 
form an operator algebra subspace of the zero-quantum operator subspace. 
This subspace is called the longitudinal magnetization and spin order oper- 
ator subspace. 

An even-order multiple-quantum operator Q e k is defined by 

Qek=Y,B kp Q2 P (p = 0,±l,±2,...) (AS) 
p 

where B kp are coefficients and the operators Q 2p are 2p-quantum operators: 

IzQ2 P \t>r) = (M r + 2P)Q 2p \^ r ) (A9) 

The definition (A8) of an even-order multiple-quantum operator shows that 
the complete set of the even-order multiple-quantum operators is a linear 
subspace of the Liouville operator space. Here will prove further that all the 
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even-order multiple-quantum operators form an operator algebra subspace 
of the Liouville operator space. First of all, the product operator of any two 
2p-quantum operators is an even-order multiple-quantum operator. It can 
be found easily from Eqs.(Al)-(A4) that 

IzQ2 q Q2 P \^r) = [M r + 2(p + q)}Q 2q Q 2p \y r ) (A10) 
This equation indicates that the product operator Qi q Qi v is a 2(p + q)- 
quantum operator, i.e., an even-order multiple-quantum operator. It turns 
out from Eq.(AlO) and the definition Eq.(A8) of an even-order multiple- 
quantum operator that the product operator Q e kQ e i of any two even-order 
multiple-quantum operators Q e k and Q e i is still an even-order multiple- quan- 
tum operator. Therefore, all the even-order multiple-quantum operators can 
form an operator algebra subspace of the Liouville operator space. 

Obviously, it follows from the definition Eq.(A8) of an even-order multiple- 
quantum operator that the even-order multiple-quantum operator subspace 
contains the whole zero-quantum operator subspace. 

There are some important properties of the even-order multiple- quantum 
operator subspace. One of which is that the exponential operator exp(±i\Q e ) 
constructed with a Hermite even-order multiple-quantum operator Q e is still 
an even-order multiple-quantum unitary operator and can be expressed as a 
sum of base operators {Q e k} of the operator subspace 

exp(±zAQ e ) =£ f p {±\)Q ep . (p = 0, ±1, ±2, ...) 
v 

where / P (±A) are coefficients. Another is that any even-order multiple- 
quantum operator is transferred into a sum of base operators of the operator 
subspace when it is acted on by an even-order multiple-quantum operator. 
These properties are obviously a direct consequence of the closed property 
of the even-order multiple-quantum operator algebra subspace. 

The properties of the longitudinal magnetization and spin order, the zero- 
quantum, and the even-order multiple-quantum operator subspace of the Li- 
ouville operator space of a two-state quantum system like a spin system may 
be helpful for simplifying the decomposition of the time-evolutional propaga- 
tor and the determination of unitary time evolution of the quantum system, 
and the decomposition of unitary transformations of qauntum computation 
into a sequence of one-qubit gates and the two-qubit diagonal quantum gates. 

Appendx B 

In a cold trapped ion system the two-qubit diagonal quantum gate 
Rmn(^mn) = ex P ( _ ^ A mri 2/ mz J n2 ) may be constructed by six laser pulses. This 
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elementary gate is really prepared with indirect interaction between a pair 
ions in the system since the interaction is set up by an intermediate media, 
i.e., phonon, while the direct interaction occurs between ions and phonons in 
the system. Adopting Cirac and Zoller s notation (Cirac & Zoller,1995), the 
elementary gate for a pair of ions m and n can be explicitly prepared by the 
following laser pulse sequence: 

Rmn{\mn) = Cft° (0 3 ) (&) K'' (fc) ^ (*l) V? (01 ) tft° (0o) 

where the parameters <^ (i = 1,2) and 6i (i = 1,2) are the laser phases 
applied to the ion n and 4> i (i = 0, 3) are the ones applied to the ion m, and 
they are not independent but subjected to the following relations: 

Q - 3 = 7T + 2(0! - 2 ), 6 1 - 9 2 = TV + 4(0! - 2 ). 

The parameter A m „ can be determined by 

A m „ = 2tt - 2(0 X - 02 ) 
Therefore, by adjusting suitably the phase difference (0 X — 2 ) of the lasers 
applied to the ion n the desired parameter X mn can be obtained. 

Figure 1. The preparation for the direct two-body elementary propaga- 
tor: SE n = Rki(Xki) = exp(-i\ M 2I kz Ii z ) (X M = 4 n nJ H T) in a coupled N-spin 
(1=1/2) system with Hamiltonian of Eq.(12). The N spins except spins k and 
I are divided into several groups {p}, {s}, {w}, where both two coupled 
spins k and I are not coupled with those spins of group {p} and there may 
be interaction among these groups but is not coupling between any arbitrary 
two spins in each of these groups. The spin echo sequence SE\ with selective 
180° RF pulses applied simultaneously to the two spins k and I and all the 
spins of group {p} refocuses their chemical shifts (one-body interactions) and 
undesired two-body interactions with any other groups but leaves the desired 
two-body interaction 2Ik z hz- The second sequence SE 2 with four SEi units 
and selective 180° pulses applied to spins of group {s} refocuses further one- 
and two-body interactions of these spins with the rest groups. Finally in 
sequence SE n all the undesired one- and two-body terms are refocused but 
only the desired term 21 kz l iz is retained and hence Rki{Xki) is obtained. 
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